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Abstract: In this paper we study dynamics of the closed FRW model with 
holonomy corrections coming from loop quantum cosmology. We consider mod- 
els with a scalar field and cosmological constant. In case of the models with 
cosmological constant and free scalar field, dynamics reduce to 2D system and 
analysis of solutions simplify. If only free scalar field is included then universe 
undergoes non-singular oscillations. For the model with cosmological constant, 
different behaviours are obtained depending on the value of A. If the value of 
A is sufficiently small, bouncing solutions with asymptotic de Sitter stages are 
obtained. However if the value of A exceeds critical value Ac = — 21mpj 

then solutions become oscillatory. Subsequently we study models with a massive 
scalar field. We find that this model possess generic inflationary attractors. In 
particular field, initially situated in the bottom of the potential, is driven up 
during the phase of quantum bounce. This subsequently leads to the phase of in- 
flation. Finally we find that, comparing with the flat case, effects of curvature do 
not change qualitatively dynamics close to the phase of bounce. Possible effects 
of inverse volume corrections are also briefly discussed. 
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1. Introduction 

In the recent years methods of the Loop Quantum Gravity (LQG) |l| have been 
successfully applied to quantise cosmological minisuperspace models. General 
prediction of this approach, called Loop Quantum Cosmology (LQG) [@], is avoid- 
ance of the initial singularity. Namely in LQG singularity is replaced by the phase 
of bounce. This result seems to be generic for the homogeneous cosmologies and 
was confirmed in the numerous investigations. Loop quantisation of the homoge- 
neous and isotropic FRW models have been performed in Ref. p, ||, |[ {K = 0), 
Ref. 0] {K = 1) and Ref. @, |^ {K = —1). In the papers cited above strict, 
fully quantum considerations were performed. However models only with a free 
scalar field and cosmological constant were studied. 

Dynamics of the models also can be traced with the effective equations of 
motion. Then classical dynamics is modified with certain nonperturbative quan- 
tum gravitational corrections. These corrections becomes unimportant for the 
low energy densities leading to the classical equations of motion. However for the 
scales comparable with the Planck energy densities corrections produce effective 
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"quantum repulsion" leading to the bounce. Studies of the effective dynamics 
have recently attracted substantial interest. It is mainly due to the fact that it 
gives the access to studying various physical consequences of the loop quantisa- 
tion. Moreover, semi-classical considerations are in many points simpler than the 
fully quantum treatment. Often approximate solution can be found before the 
strict quantum results are obtained. In particular models with self-interacting 
scalar field can be easily studied in this effective formulation. Therefore both 
fully quantum and effective methods are complementary to each other. Since 
now effective flat FRW [|1^, |11| and Bianchi I |[T^ models in LQC were studied. 



In the present paper we study effective dynamics of the closed {K = 1) FRW 
model. 

We have to stress that alternative effective description is developed by Bo- 
jowald and his collaborators. In the Bojowald's approach information about dis- 
persions and higher moments of the observables can also be obtained. Recently 
this approach has been applied to Bianchi I model |T^. However we do not apply 
this approach in our considerations. Recently also new approach to quantise loop 



quantum cosmological models has been proposed [|14| , p!5[| . This formulation leads 
to the bouncing solutions but no effects of so called Thiemann's term, which 
leads to inverse volume corrections, are present. Obtained effective dynamics is 



therefore the same like this studied in |T^, |T3] where only holonomy corrections 
were taken into account. 

Canonical variables for the closed FRW model are similar like in the fiat case. 
Namely original Ashtekar variables {A\, E^) are parametrised by the conjugated 
(c, p) variables with the Poisson bracket 

{c,p} = ^—- (l-l) 

Parameters (c, p) can be expressed in terms of the standard FRW variables in the 
following way 

p = ( —] and c = 7a— (1.2) 
\aoJ ao 

where ao = 2, fiducial volume Vq = 27r^ao = 167r^ and Iq = Vq^^ = (167r^)^/^ = 
(27r^)^/^ao. It will become useful when correspondence with classical case will be 
studied. 

In this paper value of Barbero-Immirzi parameter 7 = 0.2375, following Ref. 



16 , is assumed. We also use notation k = SnG. 
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2. Effective equations for the closed FRW cosmology 



Starting point of our considerations is the effective Hamiltonian 



n 



cflf 



sin I /ic — /i— - — sin I fi— + (1 + 7 



lo 



4 



+ Hm (2.1) 



which was derived in Ref. 0. Here /i is quantisation parameter and its form 
depends on the particular quantisation scheme. In this paper we choose so called 
/i-scheme, then 



(2.2) 



where A = 2v^7r7/pj. Earlier, /ig-scheme was also used, however it was realised 
later that it has wrong classical limit. Recently other quantisation scheme was 
proposed [l^] and studied in Ref. [|l^, |l9l . In this scheme new types of finite scale 



factor singularities appears. Another possibility is those presented in [0, |T5 
where Ji = X/-^/p and A is some unknown constant unrelated with A. 

As a matter content we consider homogeneous scalar field with the Hamilto- 



nian 



(2,3) 



With this assumption it will be possible to study realisation of the inflationary 
phase in the bouncing universe. Energy density of the scalar field can be derived 
from the following expression 



P 



1 dHrn 

p3/2 Qp^ ■ 



(2.4) 



Since the Hamiltonian constraint = is fulfilled we obtain 



2 1- -lo 
sm [fic-fi- 



2 I -lo 
sm [fi- 



-(1+7 



2\h'' ''0 



+ 



(2.5) 



where we have defined critical density 

3 VS 



Pc 



In the rest of the paper we choose the gauge = 1. 



0.82ppi. 



(2.6) 
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The equations of motion can be derived with the use of Hamilton equation 



/ = {/,Kff} 

where the Poisson bracket is defined as follows 



{f,9} 



d£dg_ _ d£dg_ 
dc dp dp dc 



+ 



df dg _ df dg 

d(j) dp^ dp^ d(j) 



From this definition we can retrieve the elementary brackets 

{c,p} = ^^^ and {0,p^} = l. 
Based on Eq. |2]^ we derive equation for the scalar field part 



P4> = W.Kff} = -P^''^'^ 



dV 



what leads to the classical equation 



■1 3Pl dV ^ 
2p dcp 



(2.7) 
(2.8) 

(2.9) 

(2.10) 
(2.11) 

(2.12) 



Here equation of motion for the scalar fields does not feel the quantum corrections. 
For the gravity part we calculate 

87rG7 9Kff 



what together with Eq. p.5| give us effective Friedmann equation 



p_ 
2p 



where 



Piip) = -P' 

P2{P) = P 



sm 



SttG 1 
3 pc 



lAlo 
p 2 



(P-P1(P))(P2(P) - P) 



(2.13) 



(2.14) 



;i+7 



A/2 
Ap 



na 



2' 



1 — sin 



'A/o 
p 2 



A/2 
2^ ^'■0 

4p 



Pc + 



(2.15) 
(2.16) 



It is clear from the above expressions that relation P2{p) — Pi{p) = Pc is fulfilled. 
Approximations performed in ( p.l5[ ) and ( |2.16| ) are vahd for p > ^ ^ 18.9 l^.^. 
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Alternatively expression for the modified Friedmann equation can be written 
as follows 

H^ = ^p^^ (2.17) 
3 a-^ 

where effective constants are expressed as follows 



GeS — G 



Pc 



Pc 



3 V'o/ Pc 



(2.18) 
(2.19) 



Expression for the Kes does not depends on the energy density p. Therefore it 
can be studied without specifying form of the matter. We expect that classical 
limit 1 for p —* oo should be recovered. It is in fact a case, what it is clear 

from Fig. |l|. For p ^ 0, factor K^.^ grows and can potentially play important 



Figure 1: Effective parameter of the curvature K^s. 

role. We study this issue in the subsequent sections. 

Defining energy density and pressure of the scalar field 



P 



(2.20) 
(2.21) 



we can derive equation 



J TJ 

— = -AttG Hp + v)- M{p)p, 



1 + 2^-2^ 

Pc Pc. 



(2.22) 
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where to simplify notation we have defined 



M{p) = {l + Y) 



A/2 



6p 



— W — sm 
3 y p 



— - I cos 

p 2 



'A/o 
p 2 



(2.23) 



It is worth to mention that taking /o = in (|2.14|) and (|2.22| ) we obtain 



StiG f p 

-irp 1 — 

3 V Pc 



2^ 



(2.24) 
(2.25) 



recovering the case K = 



3. Models with A and free scalar field as 2D dynamical 
systems 

Effective loop dynamics of the flat FRW models with A and free scalar fleld 
was studied in |TT|]. In this reference, fully analytical solutions of the model 
were found. Unfortunately in case of the closed model the equations become 
transcendental and cannot be solved analytically. Therefore other approach, to 
investigate dynamics, have to be applied. In our considerations we use methods 
of the qualitative analysis of the dynamics. Useful technique applied here is 
compactification of the phase space. Based on this approach it will be possible 
to study global properties of dynamics. 

In this section we consider models with cosmological constant and free scalar 
field. For these models dynamics reduce to 2D dynamical systems, what simplify 
qualitative analysis. Moreover one has to remember that physical solutions are on 
the surface of Hamiltonian constraint Ti.es{p, p) ~ 0. Therefore domain allowed for 
the motion reduce to one dimensional subspace. It means that on 2D phase plane 
{p,p), only one trajectory represents physical motion. Other possible trajectories 
does not fulfil 7ies{p,p) ~ 0. It is also interesting to note that in all cases, the 
same maximal value of Hubble parameter -ffmax is reached. In fact this is the 
same -ffmax = a/^Pc as this obtained for the fiat models. 

3.1 Cosmological constant 

In the first case we consider model with cosmological constant only. Energy 
density and pressure of the cosmological constant are given by 

A A 



Dynamics of the system is governed now only by modified Friedmann equation 
( p.l4| ). Since left side of this equation is > on can find physical domain of 
motion form this condition. We plot function given by Eq. ( p.l4|) in Fig. ^ 
We find that while pa < Pc then domain for the physical motion is bounded from 




two sides. Therefore we expect oscillatory solutions in this regime. While pa = Pc 
the upper bound reach infinity. In turn, when pA > Pc upper bound disappears. 
It is useful to define critical value of cosmological constant 



A. 



PI 



27r7^ 



21m 



PI 



(3.2) 



since bifurcation of solutions occur when value of cosmological constant cross Ac. 
This critical value was introduced in Ref. [H] while studying flat models in LQC 
It was show there that there is no physical solutions for A > A^.. In the curved 
models we however expect physical motion in this region. 

Now we can investigate dynamics of these three cases on the phase phase 
portraits shown in Fig. |[ We see tree mentioned types of solutions. It is worth 
to stress that only thick lines represent physical trajectories. Other trajectories 
shown in Fig. ^ are not on the surface of Hamiltonian constraint 'Hcs{p-,p) ~ 0. 



Based on approximations performed in (|2.15|) and (|2.16|) we find 



A' 



3 ll 
4A- A, 



A- Ac 



(3.3) 
(3.4) 
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Figure 3: Phase space diagrams for the model witli the cosmological constant. Shaded 
regions are unphysical because of < 0. Case a) < pA < Pc, b) pA = Pc, c) 
PA > Pc- Physical solutions, on the surface of Hamiltonian constraint TCesiPjP) ~ 0; 
are distinguished by thick (red, blue, green) lines. 



Expression for Pmax is however valid only for A > A^,. When A < Ac we have 
Pmax oo. In this case asymptotic de Sitter stage is obtained as in the flat case. 

3.2 Free scalar field 

Now we get on to example of free scalar filed which is widely used in LQC It 
is due to the fact that, in this model, scalar field can be treated as internal 
clock. It is of great importance when purely quantum cosmological models are 
constructed. However it is not necessary on the semi-classical level. Therefore 
models with no well defined global internal time can also be studied. Example 
of this will be given in the subsequent section when model with self-interacting 
field is studied. However even in that case well defined intrinsic time can be 
introduced at intervals. 

Considered model with free scalar filed was studied before in Ref. [^]. In 
this reference fully quantum analysis of the model has been done. In particular 
wave function for this model was computed. It was also shown that semi-classical 
dynamics follows precisely the quantum evolution. It is important result approv- 
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ing applicability of the semi-classical approach. Here we would like to complete 
earlier investigations. 

It is worth to remind here that energy density and pressure of the free scalar 
field are given by 

P^ = ^ and n = (3.5) 



where is constant of integration here. Based, as earlier, on Eq. ( p.l4| ) we can 
find the regions > 0. This gives us regions allowed for motion. We show 
exemplary curves in Fig. ^. It can be seen that doubly restrict region is allowed 
for the motion. This region is enlarged with increasing value of parameter p^. 




Based on analysis of positivity of function one expect that solutions are 
represented by oscillations. Namely universe bouncing between two edges of the 
region allowed for motion. It can be approved looking on phase portrait for 
this system, shown in Fig. |^. In fact physical solution is represented by closed 
curve. Similarly like in the previous case, physical solution lie on the surface of 
Hamiltonian constraint. Here we can also find 



n2 




Pmin ~ 1/ ^ and pmax ~ \l ^riPl (3-6) 



are vahd for p ^ ^ 18.9 Zpj. 



based on approximations performed in ( p.l5[ ) and ( p.l6[ ). These approximations 
A/2 

4 
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Figure 5: Phase space diagram for tlie model with free scalar field. The shaded regions 
are unphysical because < there. Physical solution, on the surface of Hamiltonian 
constraint 7ies{p,p) ~ 0, is distinguished by thick (red) line. 



4. Qualitative analysis of the model with quadratic poten- 
tial function 

Self-interacting scalar fields are commonly used in the modern cosmology. It 
comes from the fact that violation of the strong energy condition p + 3p > can 
be obtained with them. Namely p<^ + 3p<^ = 2(0^ — V{(f))) and it is possible to 
have V{(j)) > 0^ what leads to p + 3p < 0. In consequence accelerating stage of 
the evolution of the universe occurs. This mechanism can be therefore applied 
to modelling dynamics of the inflation or present era of dark energy. Here we 
concentrate on this former. It is due to the fact that stage of inflation is proceeded 
by the Planck era and both phases should be dynamically related. Namely initial 
conditions for the inflation should be given in the Planck epoch. In our approach 
the Planck epoch is described by the closed loop cosmology. It is interesting to 
see how phase of inflation emerge in this scenario. In fact this was already studied 



in case of the flat model |T^. It was shown that phase of bounce can be naturally 
connected with the inflation. Moreover inflationary attractor is generic and is 
realised for the broad range of the initial conditions. Here we are going to check 
whether curvature can affect this scenario. 
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In our considerations we assume that considered scalar field is massive. Then 
energy density is given by 



(4.1) 



Dynamical system, in case of self-interacting scalar fields, can be written in the 
form 



= 2/, 

3 z dVi(p) 

2p dcp 

p = z, 



(4.2) 
(4.3) 
(4.4) 

z = 2^p.p|2(sin^ (^|) - (1 + f )^ + 1(1,^ + T/(0))) X (4.5) 

_ A/2 1 .1 _ . N 

(4.6) 
(4.7) 
(4.8) 

■^P Pc ^ /J 

together with the constrain 



p 2' ' ■ ' ' 4p ' pc 

, .(i. 

p 2' ' ' ' Ap pc^2' 



1/1 /A, . , /A, , ,^ ^.Ml 3 2\ 
2V2\/]9 Up 2p pc 



X (cos(,/^/o) + (l + f)^-^(^l/^ + m) 



4p p, 
(4.9) 

The condition of the positivity of this expression implies that for a constant value 
of p the motion of the dynamical system in the phase plane (0, <p) is restricted to 



/ /1 , 2\^^0 I 2 / / ^ ^0 



))>iy2+imV>Pc((l- 



"7" )— ^— sm 



2\ ^''0 • 2 / / ^ "-O 

^ "m ( 



A/, 

Ap ^ W p 2 

(4.10) 

In Fig. ^ this region is represented by the shaded field (in this example p = 10). 
For a given time region allowed for the motion is in the form of a ring. External 
edge is due to quantum gravity effects while internal one coming from curvature. 
Radius of internal edge decreases for increasing p. In Fig. ^ we show exemplary 
trajectories on the (0, <p) phase portrait. Two general types of trajectories are 
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shown there. First type of trajectories is initiahsed in the point of maximal 
displacement of the scalar field in the potential, then = 0. Then filed falls 
down through the intermediate slow-roll regime, leading to the phase of inflation. 
This phase of inflation is seen as nearly horizontal attractor line in Fig. p. It can 
be seen also trajectories starting from lower displacements are directed towards 
this inflationary attractor. These trajectories are shown also in Fig. |^ where 
corresponding evolution on the (0, p) plane is presented. 




-2: -10 12 

Figure 6: Projection of the 3D phase space diagram for the model with the quadratic 
potential function and = 1 on the (0, (p) plane. The shaded field denotes region 
admissible for motion for constant value of p = 10, size of this region depends on the 
value of p. The red and blue dots represent location of the initial conditions taken for 
numerical simulations. 

Trajectories of the second kind are initialised in the bottom of the potential, 
where = 0. These trajectories are also directed to the same inflationary attrac- 
tor. It is worth to analyse this behaviour in more details. Considered evolution 
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is represented by the curves starting at = in Fig. |^. Here initial state is 




-2 -1 a 1 2 



Figure 7: Projection of the 3D phase space diagram for the model with the quadratic 
potential function and = 1 on the (</>, log (p)) plane. The arrows on the trajectories 
are omitted due to symmetries (every line can represent both expanding or contracting 
phase of evolution). 

given in the contracting pre-bounce universe. Field is initially in the bottom of 
the potential well. In this stage, value of the Hubble factor is neglected with 
respect to m^. Therefore field undergoes certain oscillations. It can be seen that 
is this stage field behaves like a dust matter. Universe behaves therefore like 
contracting, matter dominated one. However when bounce is approached then H 
grows and becomes significant for the dynamics of 4>. Since H is negative act as 
anti-friction term in equation of motion 

4> + 3H^ + m'^(f) = (4.11) 

driving up the potential well. For the considered trajectory, field is driven up 
to ~ ±1 in Planck units. Then typical slow-roll inflation takes place leading 
to exponential growth of p. Field, damped by the friction term, falls down the 
potential well and starts to oscillate finally. Then again matter-like stage occur. 
However in this point phase of reheating should stars. This phase however can- 
not be studied without extension of this model. Therefore physical analysis of 
dynamics should be stopped at this point. 
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5. Remarks on inverse volume corrections 



In the previous sections we have considered models with quantum holonomy 
corrections. Here we are making some remarks on the possible modifications 
due to the so called inverse volume corrections. In the fiat models definition 
of these corrections is ambiguous since it deepens on the fiducial cell volume. 
However in the considered closed model, fiducial cell volume is fixed and inverse 
volume corrections are well defined. Since effects of inverse volume corrections 
were extensively studied in literature we are not going to perform full analysis 
of them. Interested reader can find more information abut them in Ref. ^. 
In particular interesting feature of the inverse volume corrections is the phase of 
super- infiation \TL\. 



Here we are going to show where effects of inverse volume corrections start to 
be important. We begin this analysis with general effective Hamiltonian derived 
in Ref. § 



■21- -^o\ . 2 f -^o\ I /I , 2\/^^^0 

sm^ HC - /i— - sm^ /i— + (1 + 7^) 



2 V 2 ' ' ' A 



8nG-f 



-3/2 „2 



Here functions A{v) and B{v) contains inverse volume corrections and are given 
by 

^^-^ = '^2\lf\-\\\^-M~\v + l\\, (5.2) 
B{v) = K\v\\\v + l\'/^-\v-l\'/^f (5.3) 



where 



For V ^ 1 functions can be expanded as follows 



K = (5.4) 



A{v) ~ —, 5.5 

K 

B(v) ~ — . (5.6) 

V 



We show both functions ( |5.2| ) and ( |5.3| ) with corresponding approximations ( p. 51) 
and ( ^.6|) in Fig. ||. Approximations (|5.5|) and ( p.6|) are precisely those applied 
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Figure 8: Left : Inverse volume correction A(v) (straight line) together with ap- 
proximation used in this paper (dashed line). Right : Inverse volume correction B{v) 
(straight line) together with approximation used in this paper (dashed line). 

in this paper, leading to initial Hamiltonian (|2.1| ). Form Fig. ^ we see that 
corrections become important below v ^ 2. Since relation between p and v 
variables is given by 



we find p{2) ~ 2.8. Therefore for p ^ 3 inverse volume corrections becomes 
significant. So we have to be aware that results obtained in this paper are valid 
for p > 3 and below this value other quantum effects should take place. However 
only small fraction of trajectories probe this region. In fact we can always choose 
sufficient initial conditions to avoid this region. For instance for the model with a 
free scalar field, sufficiently high value of has to be chosen. Namely condition 
for p > 3 gives p^ > 18 Ipi. 

6. Summary 

In the classical cosmology positive curvature introduces —l/a^ term in the Fried- 
mann equation. In the early universe other terms like radiation l/a'* dominate 
and dynamical effect of curvature is negligible. However also quantum gravita- 
tional effects start to be important then and dynamics should be modified. Loop 
quantum cosmology gives us opportunity to study these effects. In particular in 
the effective formulation of LQC dynamics is governed by the classical equations 
with proper quantum corrections. These corrections modify also the curvature 
term — 1/a^ and therefore dynamical significance of this term can be changed 
in the quantum regime. In this paper we have investigated whether curvature 
term could play significant role in early universe due to the quantum effects. We 




(5.7) 
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have found that curvature does not affect evolution in the quantum regime. In 
particular, maximal value of the Hubble factor is the same like in the flat case, 
namely if^ax = a/^IPc- 

In our considerations we concentrated on the models with cosmological con- 
stant, free scalar field and massive scalar field. In the models with cosmological 
constant and free scalar field, dynamics reduce to 2D system. For the free field 
case, universe undergoes non-singular oscillations. This behaviour does not de- 
pend on the value of parameter p^. However for the model with cosmological 
constant only, type of solution depends on the value of A. For A < Ac bouncing 
solutions with asymptotic dc Sitter stages arc obtained. While A > Ac we find 
non-singular oscillations. Case A = Ac is intermediate state when upper boundary 
reaches infinity. Here Ac is critical value of A and is given by Ac = — 21mpi. 

Emergence of this value is due to quantum modifications of the dynamics. There- 
fore observed change of the type of solution is of the quantum origin and does 
not appear in the classical case. 

Physically interesting consequences occur in the models with self- interacting 
scalar field. Namely, the point is that, bounce gives mechanism to set good initial 
conditions for inflation. During the bounce fleld is driven up the potential well 
what is necessary to start phase of inflation. In the classical theory of inflation 
initial conditions have to be fixed suitably. Here they can be obtained from the 
different pre-bounce evolutions without need of any precise adjustment. This is 
nice feature of the bounce stage and is now confirmed for the fiat and closed 
models in LQC. 

References 

[1] A. Ashtekar and J. Lewandowski, Class. Quant. Grav. 21 (2004) R53 
[arXiv:gr-qc/0404018]. 

[2] M. Bojowald, Living Rev. Rel. 11 (2008) 4. 

[3] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. Lett. 96 (2006) 141301 
[arXiv:gr-qc/0602086]. 

[4] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. D 73 (2006) 124038 
[arXiv:gr-qc/0604013]. 

[5] A. Ashtekar, T. Pawlowski and P. Singh, Phys. Rev. D 74 (2006) 084003 
[arXiv:gr-qc/0607039]. 



- 16 - 



[6] A. Ashtckar, T. Pawlowski, P. Singh and K. Vandersloot, Phys. Rev. D 75 

(2007) 024035 [arXiv:gr-qc/0612104]. 

[7] L. Szulc, W. Kaminski and J. Lewandowski, Class. Quant. Grav. 24 (2007) 2621 
[arXiv:gr-qc/0612101]. 

[8] K. Vandersloot, Phys. Rev. D 75 (2007) 023523 [arXiv:gr-qc/0612070]. 

[9] L. Szulc, Class. Quant. Grav. 24 (2007) 6191 [arXiv:0707.1816 [gr-qc]]. 

[10] P. Singh, K. Vandersloot and G. V. Vereshchagin, Phys. Rev. D 74 (2006) 
043510 [arXiv:gr-qc/0606032]. 

[11] J. Mielczarek, T. Stachowiak and M. Szydlowski, Phys. Rev. D 77 (2008) 123506 
[arXiv:0801.0502 [gr-qc]]. 

[12] D. W. Chiou, Phys. Rev. D 76 (2007) 124037 [arXiv:0710.0416 [gr-qc]]. 

[13] D. W. Chiou, arXiv:0812.0921 [gr-qc]. 

[14] P. Malkiewicz and W. Piechocki, arXiv:0903.4352 [gr-qc]. 

[15] P. Dzierzak, J. Jezierski, P. Malkiewicz and W. Piechocki, arXiv:0810.3172 
[gr-qc]. 

[16] K. A. Meissner, Class. Quant. Grav. 21 (2004) 5245 [arXiv:gr-qc/0407052]. 

[17] J. Mielczarek and M. Szydlowski, Phys. Rev. D 77 (2008) 124008 
[arXiv:0801.1073 [gr-qc]]. 

[18] O. Hrycyna, J. Mielczarek and M. Szydlowski, arXiv:0804.2778 [gr-qc]. 

[19] L. F. Li and J. Y. Zhu, arXiv:0812.3552 [gr-qc]. 

[20] K. Vandersloot, Phys. Rev. D 71 (2005) 103506 [arXiv:gr-qc/0502082]. 

[21] E. J. Copeland, D. J. Mulryne, N. J. Nunes and M. Shaeri, Phys. Rev. D 77 

(2008) 023510 [arXiv:0708.1261 [gr-qc]]. 



-17- 



